Anomalous behavior of trapping on a fractal scale-free network 
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Abstract. - It is known that the heterogeneity of scale- free networks helps enhancing the effi- 
ciency of trapping processes performed on them. In this paper, we show that transport efficiency 
is much lower in a fractal scale-free network than in non-fractal networks. To this end, we ex- 
amine a simple random walk with a fixed trap at a given position on a fractal scale-free network. 
We calculate analytically the mean first-passage time (MFPT) as a measure of the efficiency for 
the trapping process, and obtain a closed-form expression for MFPT, which agrees with direct 
numerical calculations. We find that, in the limit of a large network order V , the MFPT (T) 
behaves superlinearly as (T) ~ I/2 with an exponent | much larger than 1, which is in sharp 
contrast to the scaling (T) ~ 1/" with 6 < 1, previously obtained for non- fractal scale- free net- 
works. Our results indicate that the degree distribution of scale-free networks is not sufficient 
to characterize trapping processes taking place on them. Since various real-world networks are 
simultaneously scale-free and fractal, our results may shed light on the understanding of trapping 
processes running on real-life systems. 



Introduction. — In the past decade, there has been a 
considerable interest in characterizing and understanding 
the structural properties of networked systems [1]. It has 
been established that scale-free behavior [2] is one of the 
most fundamental concepts for a basic understanding of 
the organization of many real- world systems in nature and 
society. This scalc-frec property has a profound effect on 
almost every aspect on dynamic processes taking place on 
networks [3], including robustness [4], percolation [5,6], 
synchronization [7], games [8], epidemic spreading [9], to 
name just a few. For instance, for a wide range of scale- free 
networks an epidemic threshold does not exist, and even 
infections with a low spreading rate will prevail over the 
entire population in these networks [9]. This is a radical 
change from the conclusions drawn from classical disease 
modeling [10]. 

In addition to the above-mentioned dynamics, some au- 
thors have focused their attention on the trapping problem 
occurring on complex networks [11-17], which is one of the 
main topic of interest for random walks (diffusion) [18, 19] . 
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The classical trapping problem first introduced in [20] is 
a random-walk issue, where a trap is located at a fixed 
position, absorbing all particles that visit it. An interest- 
ing quantity closely related to the trapping problem is the 
mean first-passage time (MFPT) , which is very important 
in the study of transport- limited reactions [21,22], and 
target search [23,24], amongst other physical problems. 
A result from previous research is that a power-law prop- 
erty can improve the efficiency of transport by diffusion 
on scale-free networks [11,15-17]: the MFPT, (T), scales 
linearly or sublinearly with the number of network nodes 
V as (T) ~ with 6* = 1 or 6* < 1, which shows that 
the efficiency of trapping processes on scale-free networks 
is even better than (at least not worse than) that on com- 
plete graphs [11], the best possible structure for a fast 
diffusion (with (T) - V). 

Although the scale-free topology has a direct effect 
on other structural characteristics (e.g., average path 
length [25]) of networks and dynamics running on them, 
it cannot reflect all the information of the network struc- 
ture. Recently, it has been discovered that many real- 
life networks, such as the WWW, metabolic networks, 
and yeast protein interaction networks have self-similar 
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properties and exhibit fractal scaling [26-28] . This fractal 
topology is often characterized through the fractal dimen- 
sion (Ib , which can be obtained by the box-counting algo- 
rithm [29,30]. It is now commonly accepted that fractal 
scaling [26] must be considered in an integral basic under- 
standing of the organization of real-life complex systems. 

As a fundamental property, topological fractality is re- 
lated to many respects of network structure and function. 
Recently, several authors have shown that the correla- 
tion between degree and betweenness centrality of nodes is 
much weaker in fractal network models in comparison with 
non-fractal models [31,32]. It has been also shown [27,33] 
that fractal scale- free networks are disassortative [34] , and 
this feature, together with fractality, makes such scale-free 
networks more robust against intentional attacks on hub 
nodes, as compared to the highly vulnerable non-fractal 
scale- free networks [27]. In addition to the distinction in 
the robustness, fractal networks exhibit lower synchroniz- 
ability than their non- fractal counterparts [33] . Although 
a lot of efforts have been devoted to fractal scale-free net- 
works [35-39] , it is still of current interest to look for a bet- 
ter understanding of the consequences of a fractal topology 
on different dynamic processes. 

In this paper, we study the trapping dynamics on a frac- 
tal scale-free network in the presence of a perfect absorber 
located at a fixed node. We obtain a rigorous solution 
for the MFPT of the unbiased random walks, which is 
computed through the recurrence relations derived from 
the network structure. The resulting formula shows that 
for a large network, the MFPT, (T), scales with the net- 
work order V as (T) ^ V^. This supcrlinear growth is 
significantly different from the linear or sublincar scaling 
previously found for nonfractal scale- free networks. 

The fractal scale-free network. In this section 
we introduce a network model defined in an iterative 
way [27], which has attracted a great amount of atten- 
tion [12,32,36,40]. We call this model iterative fractal 
scale-free tree (IFSFT). We study the IFSFT because of 
its intrinsic interest and its relevance to real-world sys- 
tems. For example, the so-called border tree motifs have 
been shown to be present, in a significant way, in real- 
life networks [41]. Moreover, the IFSFT is deterministic, 
which allows us to study analytically its topological prop- 
erties and dynamical processes taking place on it. It is 
thus a good test-bed and an ideal substrate network. 
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Fig. 1: (Color online) Iterative construction method of the 
network. The next generation is obtained by performing the 
following operation: for each edge, we replace it by a cluster 
on the right-hand side of the arrow, where each □ stands for a 
new external node, while ■ represents an internal node. 



The IFSFT, denoted by F„ after n (n > 0) iterations 
(the number of iterations is also called generation hcre- 




n=3 

Fig. 2: (Color online) The first three iterations of the network. 

after), is constructed as follows. For n = 0, Fg is an edge 
connecting two nodes. For n > 1, F„ is obtained from 
F„_i: for each edges in F„_i, two new nodes (called ex- 
ternal nodes with degree of 1) are firstly introduced and 
linked respectively to both ends of the edge; then, the 
edge is broken, another new node (referred to as an in- 
ternal node) is positioned in its middle and connected to 
both ends (sec Fig. 1). In other words, F„ is obtained from 
IFn-i by performing the following operations on every edge 
in F„_i: replace the edge by a path of 2 links long, with 
the two endpoints of the path being the same endpoints of 
the original edge, then attach a new node to each cndpoint 
of the path. Figure 2 shows the construction process for 
the first three iterative processes. 

According to the network construction, one can see that 
at each generation rii (rii > 1) the number of newly intro- 
duced nodes is L{ni) — 3x 4"i~i. From this result, we can 
easily compute the network order (i.e., the total number 
of nodes) Vn at generation n: 



Vn^J2 ^("*) = 4" + 1. 



(1) 



To facilitate the description in what follows, we distin- 
guish different nodes of F„ by labeling them as follows. For 
Fi, the newly created internal node is labeled 1, the initial 
two nodes belonging to Fq are labeled as 2 and 3, and the 
two new external nodes have labels 4 and 5, see Fig. 3. 
For each new iteration n > 1, we label consecutively the 
new nodes generated at this iteration, while we keep the 
labels of the old nodes unchanged. Namely, new nodes are 
labeled sequentially as Vn-i + 1, Vn-i -I- 2, ... , Vn- In this 
way, we label each node by a unique integer: at generation 
n all nodes arc labeled from 1 to = 4" -|- 1. 

Let ki{n) be the degree of a node i at generation n, 
which entered the network at generation rii [n-i > 0). If 
node i was an external node when it was introduced, 



ki{n) 



(2) 



otherwise, if i was an internal node when it entered the 
network, 

fc,(n) = 2"-"-+i. (3) 

From Eqs. (2) and (3), one can easily see that after each 
new iteration the degree of a node doubles, i.e.. 



ki{n) = 2 ki{n — 1) , 



(4) 
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Fig. 3; labels of all nodes of F2. 

The IFSFT presents some interesting topological char- 
acteristics. It has a power law degree distribution with ex- 
ponent 7 = 3 [32,36]. Its average path length, defined as 
the mean of shortest distances between all pairs of nodes, 
grows as a square power of the network order [32]. The 
betweenness distribution of its nodes exhibits a power-law 
behavior with exponent 7& = | [32]. Particulary, it is 
fractal with the fractal dimension c?b = 2. 

After introducing the IFSFT, in what follows we will 
study the MFPT for random walks with a single immo- 
bile trap on the network. We will show that the fractal 
property has an essential influence on the MFPT, i.e., the 
fractality can induce a general slowing down of diffusion. 

Formulation of trapping problem on the IFSFT. 

— Here wc formulate the trapping problem of a simple 
random walk of a particle on the IFSFT F„ in the presence 
of an absorbing trap positioned at the central hub node 1, 
represented as ir- To this end, we first represent F„ by 
its adjacency matrix A„ of order Vn x Vn- The entry Uij 
of A„ is cither 1 or 0: a^- = 1 if i and j are adjacent and 
Qij = otherwise. The diagonal degree matrix D„ of F„ 
is D„ = diag(fci (n), fc2(n), fci(n), A:y„(n)). Then, 
the normalized Laplacian matrix of F„ is given by L„ = 
I„ — A„, where I„ is the Vn x Vn identity matrix. 

In the trapping problem, at each time step, a particle, 
starting from any node except the trap it, moves from 
its current location to any of its nearest neighbors with 
equal probabilities. It is easy to see that in the end the 
particle will be necessarily absorbed by the trap, regardless 
of its starting location [11]. We are interested in the mean 
transmit time (first-passage time, or trapping time) T^^"^ 
for a particle, originating at node i, to first reach the trap 
ir in F„. 

Such a random walk can be described by a Markov 
chain [42] , whose fundamental matrix is the inverse of ma- 
trix B,i that is defined as a sub-matrix of the normalized 
Laplacian matrix L„ obtained by deleting from it the first 
row and column, corresponding to the absorbing node. 
The entry of the fimdamcntal matrix (B„)^^ ex- 

presses the mean residence time, which is the mean num- 
ber of visitations of node j by the particle, starting from 
node i, before trapping occurs. Thus, wc have 

Tt^=f:ibn%. (5) 

Then, the MFPT, (T)„, which is the mean of 7;^"^ over aU 
nodes distributed uniformly over nodes in F„ other than 
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Fig. 4: Mean first-passage time {T)n versus network order Vn 
on a log-log scale. The numerical results are obtained by direct 
calculation from Eq. (6), while the analytical results are from 
Eq. (25). The solid line is the guide to the eye. 

the trap, is given by 

I'm ^ „ -L . 

1—2 1—2 J — 2 

Equation (6) shows that the problem of finding (T)„ is 
reduced to calculating the sum of all entries of the funda- 
mental matrix (B„)"^. Although the expression of Eq. (6) 
seems compact, since the order, y„ — 1, of (B„)^^ increases 
exponentially with n, for large n, it becomes impossible to 
get {T)n through direct calculation from Eq. (6) as we are 
restricted by time and computer memory, and one can cal- 
culate directly the MFPT only for the first iterations, see 
Fig. 4. However, the particular construction of the IFSFT 
and the special choice of the trap location allow to calcu- 
late analytically MFPT to obtain a closed-form formula. 
The derivation details of which will be provided in the 
following section. 

Explicit expression for MFPT. — Before giving a 
general formula for MFPT, (T)„, we first establish the 
dependence of T^^'^ on n. 

Evolution scaling for trapping time. For each n, 
the values of Tj'"-* can be obtained straightforwardly via 
Eq. (5). Table 1 lists the numerical values of T!-'"'' for nodes 
for the first several generations up to 71 = 6. The numer- 
ical values quoted in table 1 show that for a given node i 
we have = 8Tj"'\ That is to say, upon growth of 

the IFSFT from generation n to n + 1, the mean time to 
first reach the trap increases by a factor of 8. This is a 
basic characteristic of random walks on the IFSFT, which 
can be established from the arguments below. 

Consider an arbitrary node i in the IFSFT F„. From 
Eq. (4), we know that at iteration n + 1, the degree of 
node i of IFSFT doubles, i.e., it grows from ki (degree 
at iteration n) to 2ki. Moreover, all these 2 hi neighbors 
arc new nodes created at iteration n + 1, among which 
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ki neighbors are external nodes, and the rest ki neighbors 
are internal nodes. We now examine the standard random 
walk in F„_|_i: Let A be the mean transmit time for a 
particle starting from node i to any of its fcj old neighbors, 
i.e., those nodes directly linked to node i at iteration n; 
and let B (resp. C) be the mean transmit time for going 
from any of the ki internal (resp. external) neighbors of i 
to one of the ki old neighbors. Then the mean transmit 
times follow the relations: 



A = i(l + C) + i(l + B), 



B = 
C = 



1 + A. 



Equation (7) has a solution yl = 8, found by eliminating 
B and C. That is to say, when the network grows from 
iteration n to iteration n + 1, the first-passage time from 
any node i to any node j (both i and j belong to F„) 
increases by a factor of 8. Thus, we have 7^*^"+^^ = 8Tj;"'\ 
which will be useful for the derivation of the exact formula 
for the MFPT below. 

Formula for the MFPT. After obtaining the scaling 
of mean trapping time for old nodes, we now derive the 
analytical rigorous expression for the MFPT. Before pro- 
ceeding further, we first introduce some notation used in 
the rest of this section. Let A„ denote the set of nodes in 
F„, and let A„ stand for the set of those nodes introduced 
at generation n. To facilitate the computation, we also 
define the following quantities for m < n: 



T 



m.tot / J 



in) 



and 



T. 



(n) 



in) 



ieAr, 



By definition, it follows that A„ = A„ U A„_i. 
we have 



(8) 



(9) 



Thus, 



rpin) _ rpin) 
n,tot n— 1 



7f (") _ O 

" )i,tot — " -t. 



^ n.tot 1 



(10) 



where the relation of T- 



(n+l) 



('1-1) 
l.tot 

t/"-* has been made use 



of. Hence, to determine T„ j.^^, one should first explicitly 
determine the quantity T^^j^j. For this purpose, we fur- 
ther separate A„ into two sets: one set of external nodes 
and the other set of internal nodes as defined in the sec- 
ond section, which are denoted as A„^oxt and A„ant, re- 
spectively. Clearly, A„ = A„^oxt U A„jnt- On the other 
hand, we have shown that the cardinality of set A„ is 
|A„| = 3 X 4^'-i, and that |A„,ext| = 2|A„ant|- Thus, we 
can obtain |A„ant| = 4""^ and |A„,oxt| = 2 x 4"~^. Then, 
two corresponding quantities can be defined: 



m.int 



E 



(«) 



T 



■(") _ \ ^ 

m,cxt / J 



T: 



in) 



(11) 



(12) 



It is obvious that 



n,tot n'mt ' n.cxt ' 



(13) 



Fig. 5: Illustration showing the relation of the mean transmit 
times for external and internal nodes. 



(7) In order to obtain T^^l^^^ one may alternatively get 
To this end, we first establish the rela- 



^ m,int ana 1 TO ext- 



tionship between the two quantities and ext- By 

construction, at a given generation, each edge connecting 
two nodes u and v (see Fig. 5) will generate three new 
nodes in the next generation: two external nodes {wi and 
W2) and one internal node (ws), and the mean transmit 
times for these three new nodes satisfy the following rela- 
tions: 



T{W2) -- 

Hence, we have 



1 + T(w), 

l + T{v), (14) 



T{wi)+T{w2)^2Tiw3). (15) 

Summing Eq. (15) over all old edges at the generation 
before growth, we can easily obtain that for all m < n, 



m.cxt m,int • 



(16) 



7fl(n) 



Equation (16) gives the relationship between r„ and 

— (ri) 

T,j , which is very significant since it is useful for the 
computation in the following text. 

Therefore, the issue of determining T„ is reduced to 

in) 

finding the quantity T„ g^t that can be obtained as follows. 
For an arbitrary external node lext in F„ , which is created 
at generation n and attached to an old node i, we have 



in) 



(17) 



since a particle starting from node icxt will be on node i 
after one jump. Note that Eq. (17) holds for any node pair 
consisting of an old node and any one of its new external 
adjacent nodes. By applying Eq. (17) to two sum (the first 
one is over a given old node and all its new external nodes, 
the other is summing the first one over all old nodes), we 
get 



T. 



|A„,ext|+ [h{n-l)xT^ 

jeA„_i 



2 

cxt ^ n — l.int 



(2t!: 



-2, cxt ^ ^ n-2.int 
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Table 1: The trapping time for a random walker starting from node i on the IFSFT for various n. Notice that owing to the 
obvious symmetry, nodes in a parenthesis are equivalent, since they have the same trapping time. All the values are calculated 
straightforwardly from Eq. (5). 



n\i 


(2,3) 


(4,5) 


(6,7) 


(8,9) 


(10,11,12,13) 


(14,15) 


(16,17) 


1 


3 


4 












2 


24 


32 


1 


13 


25 


29 


33 


3 


192 


256 


8 


104 


200 


232 


264 


4 


1536 


2048 


64 


832 


1600 


1856 


2112 


5 


12288 


16384 


512 


6656 


12800 


14848 


16896 


6 


98304 


131072 


4096 


53248 


102400 


118784 


135168 



+ (2"-^Til, + 2«-Ti:;L) 

= 2 X 4-1 + 2Ti'l,,,t + 4r!::^2,ext + • • ■ 
+2--'t[%,, (18) 

where Eqs. (13) and (16) were used. 

Now, we can determine T„ through a recurrence re- 
lation, which can be obtained easily. From Eq. (18), it is 

not difficult to write out T„_|^ as 

n+l,cxt — ^ ^ ^ "T^-'n.cxt "r ^ n-l,oxt ' +^ ^ l.cxt • 

(19) 

Equation (19) minus Eq. (18) times 16 and making use 
of the relation 7".^"+^^ = 8T-^"\ one gets the following 
recurrence relation 

ritiit=32ri"Lt-6x4r (20) 

(2) 

Considering the initial condition T2 ^^t = 168, this recur- 
rence relation is solved to obtain 

Tn;ext = ^(4"+^+3x32"). (21) 

Inserting Eq. (21) into Eq. (13) and considering the rela- 
tion r„"g^t = 2r„"„t, we have 

= I T^nL = ^ (4"+^ +3x32"). (22) 
Substituting the last expression into Eq. (10) yields 

Ti:L = 8 T^rAl + ^ (4-+1 + 3 X 32") . (23) 

Using r| j^^t = 14, Eq. (23) is solved inductively 

t(") = A (32" _ 4") ^ X 8" . (24) 
n,tot 28 ^ 56 ^ ' 

Inserting Eq. (24) into Eq. (6), wc obtain the rigorous 
expression for the MFPT for the trapping problem on the 
n-th generation of the IFSFT: 

(T^)n - 7^(8" -1) + 77^x2". (25) 



We have checked this exact solution for the MFPT against 
numerical values given by Eq. (6), sec Fig. 4. For all 
1 < n < 6, the analytical values obtained from Eq. (25) 
arc perfectly consistent with the numerical results. This 
agreement is an independent test of our theoretical for- 
mula. 

We show next how to represent MFPT as a function of 
the network order, with the aim of obtaining the relation 
between these two quantities. Recalling Eq. (1), we have 
4" = Vn - 1 and n = log4 (K - l)- Hence, Eq. (25) can 
be recast as 

(T)„ = A _ 1)1 _ 1) + I X (K _ l)i . (26) 
For a large network, i.e., Vn — > 00, 

{TU^iVn)K (27) 

with the exponent | = 1.5 much larger than 1. Thus, 
in the limit of large network order Vn, the MFPT grows 
supcrlinearly with the number of network nodes. 

Recently, it has been shown that for non-fractal scale- 
free networks with a large network order V, their MFPT 
(T) behaves linearly or sublincarly with V as (T) ~ 
with 9 < 1 [11,15-17]. However, we have seen that for the 
IFSFN, the MFPT, (r)„, increases superhnearly with F„ 
(i.e., {T)n ~ (Vn)^) irrespective of its scale-free property, 
presenting an obvious difference from the results previ- 
ously obtained for its non-fractal scale-free counterparts. 

Why is the MFPT for the IFSFT far larger than that for 
nonfractal scale-free networks such as the Apollonian net- 
work [16]? The reasons behind this discrepancy may be ex- 
plained as follows. For the Apollonian network, the large- 
degree nodes, including the trap node, are directly con- 
nected to one another and compose a core group sharing 
more neighbors, which make the Apollonian network be 
a very compact system. So, these large-degree nodes can 
be easily visited by a particle in spite of its starting loca- 
tion. The interconnection within nodes with large degrees 
makes the particle spend a short time to find the trap. 
On the contrary, in the IFSFT, the large-degree nodes are 
not linked to each other, they are exclusively connected 
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to low-degree nodes, as a resuh from the fr'actal prop- 
erty of IFSFT. In other words, there is an effective 'repul- 
sion' between the large-degree nodes in the IFSFT [27,32], 
which seems to be a main feature that distinguishes the 
IFSFT from the Apollonian network. Hence, for the trap- 
ping problem in the IFSFT, because of the isolation of 
the largc-dcgrcc nodes from each other, the particle will 
first reach a hub node, and before being absorbed it will 
spend a lot of time in the intermediate regions constituted 
by small-degree nodes, which connect indirectly the large- 
degree nodes to one another. Therefore, it takes a longer 
time for the particle to arrive at the trap. 

Conclusions. — We have studied a classical trapping 
problem performed on a dctcrministically growing scale- 
free network with fractal topology. The self-similarity of 
the network allows us to derive the exact expression for 
the MFPT. We have shown that, for large networks, the 
MFPT increases as a power-law function of the network or- 
der, with an exponent much larger than 1, in contrast with 
previous results found for non-fractal scale-free networks. 
We see that this slow transport efficiency lies with the in- 
herent fractality and its associated disassortativity. Thus, 
we can infer that scale- free networks do not always tend to 
accelerate the diffusion processes occurring on them. Our 
research may be helpful for a better understanding of the 
role that network structure plays in a trapping process. 
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